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FOREWORD 

This document presents the results of work performed 
by Lockheed-Huntsville Research & Engineering Center under 
Contract NAS8-31750 for the National Aeronautics and Space 
Administration, George C. Marshall Space Flight Center. 

The NASA-MSFC technical monitor for this study is 
John E. Key, EP42, Structures and Propulsion Laboratory. 
The performance period for this study was from 1 December 
1975 through 30 November 1976. 

Lockheed-Palo Altc Research Laboratory contributors 
to this effort were B.O. Almroth and F, A. Brogan. 
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SUMMARY 

This report contains the results of a study to develop modeling tech- 
niques to refine the STAGS SRB mathematical models. The objective was to 
increase the computational efficiency while obtaining adequate accuracy in 
the STAGS results. The parameter studies reported in Section 2 led to the 
formulation of an SRB mathematical model that is highly accurate and yet 
results in a significant savings in computer run time compared to previously 
used models. With the use of proper modeling techniques, the STAGS com- 
puter program is a valuable and efficient tool for the SRB structural analysis. 

The constructed SRB math model was used for the analysis of the 85/ 
32/5 (VV/VH/0) cavity collapse load condition for the aft end of the SRB struc- 
ture. Linear, nonlinear collapse and bifurcation buckling analyses were per- 
formed. The linear and bifurcation buckling analyses showed that structural 
integrity and positive margins of safety will be maintained for the applied 
maximum load case using a design factor of 1.25. The nonlinear collapse 
analysis predicts stresses in the skirt above the yield for a load factor of 
1.0. Applying the 1.25 design factor gives a negative margin of safety against 
the material ultimate. The predicted collapse occurs at a load factor higher 
than 1.25. 
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1 . IN TR OD UC TION 

Structural analysis of the solid rocket booster (SRB) is one of the prime 
responsibilities of NASA-MSFC's Engineering Analysis Division of the Struc- 
tures and Propulsion Laboratory. Many different finite element and finite 
difference computer programs are used by MSFC in carrying out their re- 
sponsibilities. Some of these computer programs were developed by Lockheed 
Missiles & Space Company. In particular Lockheed programs STAGS A, 
STAGS B, BOSOR, and SPAR are used extensively by'MSFC, partly because 
MSFC participated in the development and helped sponsor these programs 
and are quite familiar with them. 

In the course of their detailed work, MSFC encountered some limita- 
tions in the output of the subject programs. Since much expensive SRB 
hardware is involved it is imperative that the structural analyses produce 
accurate results. The purpose of this effort was to ensure that the subject 
programs produced accurate results for all required analysis conditions and 
applications. The latest version of the STAGS finite difference program, 
STAGS C, was installed, checked out and an SRB model formulated for the 
accurate analysis of the cavity collapse load conditions. Some discrepancies 
in the STAGS analyses performed by NASA support personnel have occurred 
in the past. These problem areas were investigated and the causes deter- 
mined and solutions obtained. 

STAGS is a computer program developed at Lockheed-Palo Alto and is 
intended for analysis of a shell type structure. The structure to be considered 
may consist of up to 30 different shell branches that are treated by use of 
finite difference discretization. A shell branch is any part of the structure 
that is thin. In addition to the shell branches, the structure may include some 
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finite elements: elastic bars, shear panels, beams and nonlinear triangular 
plate elements. 

The code can be used for: 

• Linear stress analysis 

• Geometrically nonlinear elastic stress analysis 

• Inelastic stress analysis, geometrically linear or nonlinear 

• Bifurcation buckling analysis with linear or geometrically 
nonlinear pre stress (elastic), and 

• Small vibration analysis with prestress based on linear 
or geometrically nonlinear analysis (elastic) 

• Transient response analysis, linear or geometrically 
nonlinear, elastic or inelastic. 

Any combination of point forces, line loads and distributed surface 
tractions can be applied. Loading by specification of displacements or thermal 
gradients (through the shell wall and over the shell surface) is also permitted. 
Any configuration of boundary conditions or other displacement constraints can 
be included in the analysis. Reference 2 contains the instructions for the use 
of the STAGS code. 

This report contains five technical sections; four of which describe the 
tasks performed and the final section lists several conclusions and recom- 
mendations of areas requiring further study efforts. 
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2. STAGS MODEL PARAMETER STUDY 

2.1 EFFECTS OF GRID SPACING ON LINEAR ANALYSIS 

Problems in connection with the use of the STAGS finite difference com- 
puter program have been experienced during the analysis of the SRB structure 
by NASA support personnel. Some of these problems were the result of im- 
proper grid configurations being used in the models. Parameter studies were 
made to determine the proper grid layout for accurate results in the linear, 
nonlinear, and bifurcation buckling analyses. Results concerning the proper 
grid spacing for the bifurcation buckling analysis are presented in Section 2.2 
and for a nonlinear analysis in Section 2.3. 

The aft portion of the SRB cylindrical body and the attached conical 
skirt, loaded by the cavity -collapse phenomena during water entry, are the 
structural components under investigation. The body is a ring stiffened 
cylinder of constant wall thickness reinforced with tee -stiffeners and clevis 
joints between body segments. The aft end of the cylinder has a spherical 
closure which supports the rocket nozzle. The motor case structure and 
components are fabricated from D6AC steel. A reinforcing tee-ring is 
located at the cylinder -skirt intersection. 

The conical skirt is a ring- stringer stiffened shell with skin thickness 
varying both lengthwise and circumferentially. Additional components, such 
as holddown posts and nozzle actuator support brackets, are integrated into 
the skirt structure. The skirt structure is fabricated from aluminum. A 
simplified view of this portion of the SRB configuration is shown in Fig. 2- 1 
(without the nozzle). In all analyses, (j> - 0 is the keel side and $ = 180 is the 
lee side for the water cavity collapse load cases. 
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The major discrepancy in the previous analyses is in the cone-cylinder 
intersection area. Brown personnel experienced discrepancies in stress 
values and deformations between two seemingly identical models; i.e., differ- 
ent results were obtained when the reinforcing tee -ring was modeled on the 
cylinder portion as compared to results obtained when it was modeled on the 
cone portion of the structure, all other parameters held constant. Presented 
in this section are data on the effect of variable -grid size on the accuracy of 
the analysis in this region. 

The simplified cone-cylinder model shown in Fig. 2-2 was used for the 
"grid study." The same length, radius, cone angle, etc., were maintained for 
the simplified SRB model, but a constant skin thickness was used. In addition, 
all components were assumed to be made of the same material. No stiffeners 
were added to the structure except for the reinforcing ring at the cone -cylinder 
intersection. Initially a 90 degree segment of the structure, with symmetry 
planes at the side boundaries, was used; this was later reduced to a much 
smaller segment after it was determined that the number of grid divisions 
in the circumference had no effect on the results for the internal pressure 
loading used in this study. The structure is modeled as two "branches"; one 
is the cylinder and the other is the cone. In STAGS each general shell com- 
ponent is modeled as a separate branch and the geometry, material properties 
and loading are specified for that branch. The boundary conditions at the junc- 
tion of these branches is specified by the user for compatibility in the structure. 

The model was first checked for symmetry and to determine if the cor- 
rect cylinder wall stresses were accurately computed by the program. 

Computer runs were then executed for the sample problem using a 
range of longitudinal and circumferential grid spacing s. The accuracy in- 
creased with a finer grid spacing in the longitudinal direction, i.e., more 
grid divisions per length. For each grid configuration the analysis was 
executed twice, once with the reinforcing ring on the cylinder and again 
with the ring on the cone branch. The aspect ratio, ratio of length to width. 
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Fig. 2-2 - Simplified SRB Model 
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of the grid was seen to have no effect on the results. Several grid sizes were 
investigated using both uniform and nonuniform grid point spacing. 

It was found that a variable grid spacing produced comparable results 
with the uniform grid spacing while significantly reducing computer run costs. 
This was achieved by allowing large grid spacing s in areas of small stress 
gradients, and small grid spacings in areas of large stress gradients. This 
method of grid layout is the most efficient arrangement for- minimizing com- 
puter run time for a desired level of accuracy. 

Table 2-1 shows the effect of longitudinal grid size on the results. The 
longitudinal grid spacing adjacent to the intersection is given in inches. The 
grid spacing was the same in the two branches, cone and cylinder. This 
spacing was held constant for at least two grid spacings each side of the inter- 
section. CYL denotes the analysis with the reinforcing ring attached to the 
cylinder branch and CONE denotes the analysis with the ring on the conical 
skirt. The tee-ring was modeled as a rigid stiffener. As the longitudinal 
grid size is reduced the radial displacement, W , and rotation, [3 , become 
more comparable for the two analyses. Note the results for the 2.0 and 0.5 
inch grids. 

To understand why "j3" (Table 2-1) changes sign between the "CYL and 

CONE" analyses, for any shell structure there exists a "boundary layer" at 

1/2 

the shell edges of the same order of size as the (radius x thickness) ' 
within which the shell theory does not give an accurate representation 
of the stress field. However, if the shell approximations are applicable 
in the major part of the structure, shell theory will give good estimates 
of these stresses and displacements outside this boundary layer. The 
edges at the intersection of the cone and cylinder branches are in this 
boundary layer and approximations and interpolation must be used to obtain 
compatibility and continuity with the rest of the structure. The finite differ- 
ence integration scheme used in STAGS also must be modified at the boundary. 
A "whole station" method of grid points and integration points is used to 
produce the greatest accuracy for the shell analysis. This scheme is shown 
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Reinforcing 
Tee Ring 


Table 2-1 
GRID STUDY 


Item 

Grid Size 
(in.) 

10.0 

5.0 

2.0 

0.5 

\, 

Branch 

Ring 

1 

Ol = 

ou 

CY L 

6840 

-274 

-586 

-468 

-467 




CONE 

6979 

148 

-462 

-442 


2 

CT 3 


CYL 

4168 

5494 

1994 

30 

-232 




CONE 

10600 

5765 

2108 

434 


3 

„ •V'f 

°4 

CYL 

9512 

-6042 

-3166 

-966 

-618 




CONE 

3358 

-5467 

-3032 

-1319 











4 

W 

c 

CYL 

0.305 

-0.00069 

-0.00148 

-0.00118 


5 


0 

•H 

+-> 

CONE 

0.0176 

-0.00038 

-0.00116 

-0.00118 


6 

P 

<0 

> 2 

CYL 

0.00153 

0.00447 

0.00524 

0.00236 


7 


(!) 

% 

CONE 

0.00084 

-0.00333 

-0.00447 

-0.00165 


8 

M 

X 

< 

CYL 

-60. 6 

-625 

-1043 

-1310 


9 



CONE 

-59.9 

-563 

-954 

-1225 



°l’ °Z are ■* 100 P stresses, psi. 

° 3 ’ °4 are axial bending stresses, psi 
W = radial deflection, positive outward (in.) 

_ 0-yy 

(3 = change in slope in the axial direction. 

M = axial local bending moment (in-lb) 

^ reinforcing tee -ring results with tee ring modeled as a separate 
branch. 
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in Fig. 2-3. Only one set of grid points is used but two sets of integration 
points, one for the bending strain energy (integration area A) and one for the 
membrane strain energy (integration area B). This method has been shown 
.to produce better results than the half-station integration scheme that was 
used in STAGS A. Modifications have to be made to the integration scheme 
for accurate analyses at these boundary grid points. A plot of the rotations 
on either side of the intersection for the 0.5 inch grid spacing analysis 
in Table 2-1 is shown in Fig. 2-4. The solid curves gives the rotations for 
the analysis with the reinforcing ring on the conical branch. The dashed 
curve represents the-analysis with the ring on the cylinder branch. The 
difference in the rotations calculated at the intersection is seen to be less 
than a fourth of a degree for modeling the ring on the cone or the cylinder 
branch. The calculated inflection point is within 0.20 inch of the intersection 
line by either analysis. Thus, one can. see why (5 (Table 2-1) has an apparent 
but insignificant discrepancy. 

The stresses at four points on the ring are also given in Table 2-1. 

The stresses at points 1 and 2 are the hoop tension values and 3 and 4 are 
the axial bending stresses. Runs were also made with the ring web modeled 
as a separate branch with the flange as a stiffener. This- allows bending of 
the web and less rotation of the flange than when the ring was modeled as a 
rigid member. The ring web deformation is plotted in. Fig. 2-5 and compared 
with the rigid ring rotation. The reduction in bending is noted in the stress 
values for points 3 and 4 in the last column of Table 2-1. The ring web axial 
stress value is the same as that shown in the column for the 0.5 inch spacing 
where the ring was a rigid member. 

The discrepancies between the CONE and CYL runs are also partially 
due to the interpolation using several consecutive node properties either side 
of a node to determine the displacements at that node. Where a large dis- 
continuity exists in the model, geometrical and/or stiffness, too coarse a 
grid size will not allow an accurate solution in this particular area. At the 
cone -cylinder intersection with the stiff reinforcing ring, the large grid size 
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Half-Station Scheme in STAGS A 


X 


Grid points for w and 
integration points 


Q Grid points for u,v 


* 


* 




^ Grid points (u,v,w) and 
integration points for 
bending energy 


Q Integration points for 
membrane energy 


Whole -Station Scheme in STAGS B and C 


Fig. 2-3 - Some Schemes for Two-Dimensional Finite Difference 
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Rotation 
in Radians 



Fig. 2-4 - Branch Rotations at the Intersection 
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inhibits the true radial displacement at that point. The effect of grid size on 
the radial displacement at the intersection is shown in Fig. 2-6. Two curves 
are shown; the solid graph shows the results with the same grid spacing on 
both sides on the intersection. The dashed curve shows the effect of main- 
taining a course grid in a branch on one side of the intersection while having 
a fine grid on the other branch. A fine grid spacing must be maintained in 
the immediate area of discontinuity for accurate results. , A coarse grid can 
then be established in the rest of the model with no appreciable effect on the 
results and considerable savings in runtime. 

2.2 EFFECTS OF GRID SPACING ON BIFURCATION BUCKLING ANALYSIS 

In a stiffened thin skin structure, buckling (local as well as general) 
can be the predominant factor governing design. Therefore, a thorough bi- 
furcation buckling/nonlinear collapse analysis is needed for the aft portion 
of the SRB under cavity collapse loads. The nonlinear collapse analysis will 
be discussed in Section 2.3. 

■To determine the effect of grid size on the magnitude of in-plane bi- 
furcation buckling loads the simplified cone -cylinder model described in 
Section 2.1 without stiffeners was investigated. Loading was due to a uni- 
, form axial ring load applied along the tangent to the cone at the aft end of 
the skirt. The branch lengths of the cone and cylinder were each divided 
into 11,21 and 41 rows on separate runs with .the same number of rows in 
each branch. Table* 2.2 and Fig. 2-. 7 show the computed eigenvalues l ' for 
each run. Note that the eigenvalue is converging as the grid spacing is 
decreased. STAGS B was used in this analysis. The numerical procedure 
used in STAGS B allows for convergence of eigenvalues from above with de- 
creasing grid spacing. Too coarse a grid spacing in the longitudinal direction 


‘The eigenvalue is the value by which the applied base loads must be multi- 
plied in order to give the bifurcation buckling load. In this report, eigenvalue, 
bifurcation buckling load factor, and critical load factor will be used inter- 
changeably. 
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Table 2-2 

BIFURCATION BUCKLING EIGENVALUES 


No', of Divisions 
In Each Branch 

ELg en value 

11 

9.44 

21 

5.60 

41 

4.83 



10 20 30 40 50 

Number of Divisions in Each Branch 


-7 - Bifurcation Buckling vs Grid Spacing Simplified 
Cone -Cylinder Model Using STAGS B 
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results in very high values for the critical load factor. It was determined 
that in order to achieve accurate bifurcation buckling load factors, the initial 
grid should be sized so that the half wavelength of the bifurcation buckling 
mode (that portion of the mode between inflection points) contains at least 
three intermediate grid points. In STAGS Cthe convergence of bifurcation 
buckling load factors is from below with decreasing grid spacing. However, 
the conclusions concerning intermediate grid points are still applicable. 

Several STAGS C bifurcation buckling analyses were performed to 
obtain results for the 85/32/5 cavity collapse pressure distribution described 
in detail in Section 4 of this report. The 180 deg thin skin shell model de- 
scribed in Section 2.3 was used in the analysis. Effectiveness of discrete 
stiffeners was determined by successively adding more stiffeners in analyses 
with the same grid. Finally, different grid patterns were set up to determine 
the effect of grid layout. 

With an unstiffened structure, the critical load factor was found to be 
0.08. With the two stiffeners included on the dome -nozzle interface and 
motor case -skirt intersection, the critical load factor increased to 0.11. 

For both the critical load factors mentioned above, buckling was predicted 
at the aft end of the skirt. However, in these analyses, none of the circum- 
ferential or meridional stiffeners in the skirt were included. To account 
for the meridional stiffeners an equivalent thickness of the skirt was calcu- 
lated. The equivalent thickness was determined so that the bending stiffness 
of the skin of equivalent thickness is the same as the stiffness of the skin 
with meridional stiffeners included. 

Modifying the previous structure with the three circumferential stiff- 
eners along the skirt, along with the equivalent thickness, the critical load 
factor was found to be A = 0.32, with buckling predicted in the cylindrical 
motor case just forward of the dome-motor case intersection approximately 
at Station 1800. 
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Finally, the model was modified to include the equivalent thickness in 
the skirt as well as ten discrete stiffeners: a stub, two tee rings, and two 
clevis joints along the motor case; the stiffener at the aft end of the dome; 
the main attach ring at the motor case-skirt intersection; and three stiffeners 
equally spaced along the skirt. This model includes all discrete stiffeners 
on this portion, of the SRB except the integral longitudinal stiffeners which 
were included in the orthotropic skin. Bifurcation buckling results for this 
model with different grid layouts are presented in Table 2-3. Note that grid 
spacings have a significant effect on the critical load factors and therefore 
on the predicted buckling locations. 

As pointed out earlier, in STAGS C the convergence of the bifurcation 
buckling load is from below with decreasing grid spacing. Convergence from 
below makes possible the occurrence of spurious solutions. Since the bifur- 
cation buckling load for the detailed model under the same loading was at a 
load factor of 2.5 (see Section 4) the results from the simplified model are 
suspect. Therefore, the detailed model described in Section 3 was chosen 
for all succeeding analyses. 


Table 2-3 

BIFURCATION BUCKLING LOAD FACTORS 


Grid in 
Motor Case 

Critical 

Load 

Factor 

Branch 

Station 

Angle 

(deg) 

17 x 13 EC,UL* 

0.94 

Motor 

Case 

1660 

75 

17 x 19 EC, UL 

1.20 

Motor 

Case 

1762 

180 

17 x 25 EC, UL 

1.03 

Skirt 

1880 

180 

17 x 19 UC, UL 

0.98 

Motor 

Case 

1663 

75 

14 x 13 EC, EL 

0.46 

Motor 

Case 

1630 

180 


‘EC -Equal spacing in circumferential direction, UL-Unequal 
spacing in longitudinal direction, etc. 
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2.3 NONLINEAR ANALYSIS 

A more reliable-method to determine impending collapse is to perform 
a nonlinear collapse analysis. In the nonlinear collapse analysis the desired 
range of loading is divided into intervals or increments. Each succeeding 
load increment includes the geometric nonlinearities (nonlinearities resulting 
due to large deformation theory) resulting from all preceding load increments. 
The loading from each increment is summed to achieve the desired level of 
loading. Incremental loadings can be mechanical or thermal in nature. 

Since the Univac 1108 C version of STAGS is relatively new, its non- 
linear analysis results were checked by use of the finite element program 
NEPSAP and the finite difference programs STAGS B and BOSOR4. NEPSAP 
is a Lockheed-developed three-dimensional finite element computer program 
for the nonlinear thermo -elastic -plastic and creep analysis of arbitrary struc- 
tures undergoing large deformations. BOSOR4 is a comprehensive computer 
program for the analysis of stress, stability, and vibration modes of segmented, 
ring -stiffened, branched shells of revolution and prismatic shells and panels. 

It performs large deflection axisymmetric stress analysis, small deflection 
non-axisymmetric stress analysis, modal vibration analysis with axisymmetric 
nonlinear prestress included, and buckling analysis with axisymmetric or non- 
symmetric prestress. 

A 1 deg meridional slice of a simplied SRB cone-cylinder model was 
used to check out the nonlinear analysis. The model was composed of one- 
half inch thick aluminum. The cylindrical section of the model was fixed at 
the forward end. For the analysis the reference surface of the shell was 
placed at the midplane of the thickness. Loading on the structure was due 
to a unif orm meridional line load applied at the aft end of the conical section. 

A schematic showing the model, boundary conditions, and the point of appli- 
cation of the line load is given in Fig. 2-8. Also given in Fig. 2-8 are load 
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Fig. 2-8 - Load-Deflection Curves for Nonlinear Analysis 
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displacement curves derived from, the computer codes NEPSAP, BOSOR, 
STAGS B and STAGS C. For the purpose of the discussion here no distinc- 
tion between STAGS B and STAGS C will be made since identifical results 
were obtained from these programs. In Fig. 2-8 the magnitude of applied 
loading is indicated by the line load factor LLF which is merely a multiplier 
of the^applied base load of 1500 Ib/in. (for instance LLF = 2.0 implies a line 
load of 3000 lb/in.). 

For the linear analysis indicated in the figure a small line load (LLF = 
0.4) was applied to the structure. This small magnitude of load was chosen 
to minimize the nonlinear effects in the corresponding nonlinear runs used 
for comparison. NEPSAP, BOSOR, and STAGS all produced identical results 
at this load level. The curve in Fig. 2-8 labeled linear analysis has been ex- 
tended to a value of LLF = 4.0 to indicate the magnitude of nonlinearities that 
exist in the other runs. 

Collapse loads indicated by STAGS were at a line load factor of approx- 
imately LLF =5.4 whereas BOSOR and NEPSAP indicate collapse in the range 
of LLF = 4.4 to LLF = 4.8. Several things affect the load level at which collap 
occurs. The type of terms included in the energy expressions and the numer- 
ical method of solving the resultant equations can affect the level at which 
collapse is indicated. However, probably the most significant reason for 
differences in the collapse level is the method by which STAGS automatically 
reduces the size of the load increment when convergence within a load incre- 
ment is difficult. For both BOSOR and NEPSAP a constant load increment 
of LLF = 0.4 was used. For impending collapse, when the load increment is 
too large, numerical instabilities can occur, causing erroneous results for 
this last load increment. Difficulty in convergence occurred for BOSOR and 
NEPSAP at LLF = 4.0, LLF = 4.4, and LLF = 4.8. For BOSOR the displace- 
ment pattern at LLF = 4.8 was completely different from previous solutions, 
probably representing a solution in the postbuckling range. A BOSOR bifur- 
cation buckling analysis was not performed to determine the corresponding 
eigenvalue (LLF) and buckling location. For NEPSAP, the displacement 
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pattern was consistent for LLF = 4.8; however, the incremental displacement 
in going from LLF = 4.4 to LLF = 4.8 was approximately ten times the total 
displacement that had occurred at the end of the preceding step, indicating 
the onset of collapse. For the STAGS nonlinear analysis, when convergence 
was difficult, the size of the load increment was automatically cut in half and 
then halved again when necessary to achieve convergence. Within the whole 
analysis using STAGS, the size of the load increment decreased from a value 
of LLF = 0.5 in'the early stages to a value of LLF = 0.031 in the latter stages 
of the analysis. By using this procedure STAGS was able to achieve conver- 
gence up to a value of LLF = 5.406 at which maximum allowable execution 
time for the analysis was reached and analysis was then terminated. Due to 
the relatively small load increment needed in STAGS to achieve convergence 
it is evident that collapse is impending. Indicated in Fig. 2-8 is a STAGS 
bifurcation buckling load of LLF = 5.6. The STAGS nonlinear analysis is 
asymptotic to this buckling load, thereby showing compatibility between 
collapse and buckling analysis. 

It should be mentioned here that at no time did there exist any discrep- 
ancies between the nonlinear STAGS results and the linear STAGS results as 
NASA support personnel experienced. For their analysis, there were signif- 
icant differences in the slopes of the load displacement curves between the 
linear analysis and the initial load increments of the nonlinear analysis. In 
order to determine the reason for discrepancies NASA support personnel 
experienced, exact models and load conditions used by them would have to 
be studied. The end results of the nonlinear analysis parameter study was 
to derive a simplified model with computer run times significantly less than 
those for the detailed model described in Section 3. To this end a simplified 
model was constructed. 

The 180 deg simplified model includes all stiffeners in the detailed 
model with the exception of the hold-down posts. Ten discrete stiffeners 
in all have been included; a stub, two tee-rings, and two clevis joints along 
the motor case; the stiffener at the aft end of the dome; the main attach 
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ring at the motor case-skirt intersection; and three stiffeners equally spaced, 
along the skirt. The following thickness and materials were used for each 
branch: 

Motor Case — Steel, Thickness = 0.5 in. 

Dome —Steel, Thickness = 0.362 in. 

Skirt —Aluminum, Thickness = 0.5 in. 

In addition, orthotropic material properties were used in the skirt to account 
for the longitudinal stiffeners. 

The 90/45/10 cavity collapse pressure distribution described in the 
initial NASA loads document was used in the analysis. The same pressure 
distributions exist for the 85/32/5 load case described in Section 3 but at a 
reduced magnitude. Undeformed plots of the simplified model are given in 
Figs. 2-9 and 2-10. These plots were obtained from the NEPSAP program. 

A nonlinear collapse analysis was performed using the aforementioned 
simplified model. The results at Station 1800 for <f) = 180 deg (this is the 
point in the structure with greatest nonlinearity in load path) are shown 
in Fig. 2-11. Also shown in the figure are the results from the detailed model 
described in Section 3. The initial slopes of the load-displacement curves for 
the two models differ for two reasons: (1) the hold-down posts and thickened 

skin portion of the skirt have not been included in the simplified model, and 
(2) a thickness of 0.5 was used in the simplified model whereas a thickness 
of 0.52 was used in the detailed model. Due to the simplifications for the 
parameter study, the simplified model is slightly more flexible than the de- 
tailed model. This is evident on inspection of Fig. 2-11. Still, displacements 
differ by less than 15% for the two models. Convergence in the detailed model 
at a load factor greater than 1.5 was not obtained, indicating the onset of col- 
lapse. However, in the simplified model, convergence did not become a 
problem until a load factor of 1.7 was reached. 
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ORTHOGRAPHIC VIEW Y-Z PLANE 



Fig. 2-9 - Orthographic View of Simplified SRB Thin Shell Model 
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TOTAL STRUCTURE 
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Fig. 2- 10 - Isometric View of Simplified SRB Thin Shell Model 
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Fig. 2-11 - SRB Motor Case, Station 1800, $ = 180 deg, 90/45/10, 
Cavity Collapse Loads 
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The s imp lified model described above offered reasonably good overall 
load -displacement characteristics. However, to determine the structural 
integrity of the SRB these overall results are not enough. Very detailed 
stresses within the skin are needed and these could not be determined using 
the simplified model. For this reason, the detailed model described in 
Section 3 was used for all succeeding analyses. 

2.4 REFERENCE SURFACE LOCATION 

In the STAGS program, both versions B and C, a thin shell is located 
by defining a reference surface and the offset from the reference surface 
to the midplane of the shell. NASA support personnel have determined dis- 
crepancies dependent upon the way in which the midplane of the shell is de- 
fined. A series of STAGS C computer runs was made to determine the effect 
of reference surface location on results. 

The simplified unstiffened cone-cylinder model described in Section 
2.3 was used in this study. Linear responses of the structure were deter- 
mined. The reference surface was located at the midplane of the shell and 
also at the outer surface of the shell (skin centerline offset at a distance of 
one-half the thickness from the reference surface). 

Under internal pressure loads identical results were obtained, regard- 
less of the location of the reference surface. For the reference surface 
located at the midplane of the shell wall and line loads applied at the aft end 
of the skirt, the calculated stress resultants along the aft end of the skirt 
were identical to the applied line loads. However, under the same loadings, 
with the reference surface located at the outer skin surface, significant 
variations in the stress resultants from the applied line loads occurred at 
the point of load application. It was determined that line loads in STAGS B 
and STAGS C are applied at the reference surface, not at the skin centerline. 
Therefore when the reference surface is not at the skin midplane, the meri- 
dional line load plus a moment to compensate for the eccentricity must be 
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applied. When applying this equivalent load system with the reference sur- 
face at the outer surface of the shell, the results achieved were identical to 
those with the reference surface at the midplane of the shell. 

It should also be pointed out that the discrepancies occurring due to 
reference surface location were quite localized and completely damped out 
within the ’’boundary layer" region defined by the meridional distance (thick- 
ness x radius)^^. 

2.5 IMPERFECTION SENSITIVITY 

Imperfections from manufacturing tolerances could possibly affect the 
load carrying capability of the aft portion of the SRB. A tolerance of +0.5 
inch on the overall diameter of the motor case is specified on the part draw- 
ings. It is apparent that slight variations in initial geometry would have 
little effect on the response of the structure for nonsymmetric loading 
patterns. However, for the symmetric pressure distributions that exist for 
vertical entry of the SRB on water impact there is a possibility that imper- 
fections could significantly affect the collapse loads of the structure. 

The STAGS C computer program handles imperfections by a user 
written subroutine, WIMP. Basically, this subroutine defines variations 
in normal displacements from the original reference surface due to imper- 
fections. The prescribed structure (original structure with imperfections) 
is assumed to be in a stress free state, and a nonlinear analysis is performed 
on this stress free structure. 

In order to study the effect of imperfections in the motor case the 
following variations in normal displacements along the circumference were 
assumed: 

Wj_ (0) = - ~ cos. 26 0 < 0 < 180 deg 
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where T is the specified tolerance on the diameter. In the motor case T — 

0.5 inch. Variations along the axis of the motor case were of the form 

X. < X < X. 

1 - - 3 

where x. and x. are axial coordinates at the beginning and end of the axial 
1 3 

imperfection. The total variation in normal displacements was then 

W (6,x) = Wj (6) W 2 (x) 

Note that at the end points for the imperfections (x = x^, x = x^, 0 = 0 deg, 

0 = 180 deg) the slopes are zero so that no discontinuities exist in the skin. 

Using the STAGS C program, the imperfections as defined above were 
applied to the detailed SRB model described in Section 3 of this report. The 
end points for the longitudinal variation in imperfections were chosen to be 
at the tee rings on the motor case. This corresponds approximately to the 
point of maximum pressure for the 90/0/0 load case described in Revision 
B to the SRB loads document (Ref. 1). Pressure distributions for the 90/0/0 
load case were applied to the structure. 

Results of the linear and nonlinear analysis with and without imperfec- 
tions are shown in Fig. 2-12. The lee side (0 = 180 deg) radial displacement 
midway between the tee rings is plotted against the load factor. The load 
factor is the multiplier for the 90/0/0 pressure distribution defined in the 
loads document (i.e., a load factor of 1.0 defines the 90/0/0 loading condition 
and a load factor of 2.0 doubles this). Comparison of the results for the 
model with and without imperfections indicate that the magnitudes of imper- 
fections in the motor case due to manufacturing tolerances do not significantly 
affect the response of the structure. In fact, at a load factor of 2.6 the dis- 
placements were within 1% of each other. 
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Fig. 2-12 - Motor Case Radial Displacement on Lee Side, Midway 
Between Tee Rings for 90/0/0 Load Case 
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The manufacturing tolerances on the SRB skirt are of the same order 
of magnitude as those in the motor case. Also, the pressure loads on the 
skirt for the 90/0/0 load case are smaller than those on the motor case. 
Since imperfection sensitivity was not detected in the motor case, it is antic- 
ipated that the skirt will exhibit the same characteristics. 

The load factor value at which nonlinear collapse occurs might be sig- 
nificantly affected by imperfections, but a load factor of a least 3.0 would be 
reached before collapse. Load factors beyond 3.0 are of little interest in 
the determination of the structural integrity, of the SRB (see Section 4). 
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3. STAGS C SRB MATHEMATICAL MODEL 

3 . 1 STRUCTURE D ESCRIPTION 

The aft portion of the SRB cylindrical body and the attached conical 
skirt, loaded by the cavity -collapse phenomena during water entry, are the 
structural components to be modeled- The body is a ring stiffened cylinder 
of constant wall thickness reinforced with rectangular stiffeners, tee- stiffeners 
and clevis joints between body segments. The aft. end of the cylinder has a 
hemispherical closure which supports the rocket nozzle. The motor case 
structure and components are fabricated from D6AC steel. A reinforcing 
tee-ring is located at the cylinder -skirt intersection. 

The conical skirt is a ring -stringer stiffened shell with skin thickness 
varying both lengthwise and circumferentially. Additional components, such 
as holddown posts and nozzle actuator support brackets, are integrated into 
the skirt structure. The skirt structure is fabricated from 2219-T87 alumi- 
num. A simplified view of this portion of the SRB configuration showing the 
major geometry dimensions is shown in Fig. 3-1 (without the nozzle). 

The data have been compiled from the SRB drawings and coded into the 
STAGS C format. The skirt data were taken from NASA-MSFC 14A00 305 
(12 January 1976) drawing and the SRM data from Thiokol IU50185 and IU50129. 

3.2 SRB MATHEMATICAL MODEL ■ 

The math model is a 180 degree symmetrical structure composed of 
four branches: Branch 1 is the motor case cylinder from the reinforcing 
ring at Station 1613.50 to the aft dome tangent point at Station 1828.85; Branch 
2 is the aft dome; Branch 3 is the remaining portion of the cylinder between 
the aft dome tangent point and the cone-cylinder intersection; Branch 4 is the 
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conical skirt. The tee-rings and clevis joints are included on Branch 1 and 
the heavy nozzle attach ring on Branch 2. The large reinforcing tee-ring at 
the cone -cylinder intersection is included on Branch 3. The three internal 
rings and the longitudinal stiffeners are included on the skirt as well as the 
holddown posts and actuator support brackets. The heavier primary stiffeners 
that tie into the rings at the aft end of the stiffeners are modeled as discrete 
stiffeners and the lighter intermediate stiffeners are smeared into the skin. 
The skin then becomes orthotropic with increased stiffness in the longitudinal 
direction. This maintains the structural stiffness in the model but reduces 
the number of grid points required in the circumferential direction. This 
model has 3944 nodes and 6587 equations. 

The variation in wall thickness in Branches 2 and 4 are input using 
Subroutine WALL. The subroutine for inputting stiffeners off the gridlines 
has not been completed and checked out for STAGS C, so the holddown post 
diagonal members had to be input point-by-point on the gridlines. The stiff- 
ness properties were prorated between the affected nodes. 

Two SRB math models are used for the cavity collapse analyses. They 
are designated as the "original" and the "rotated" model. The "original" SRB 
math model used for the analyses is oriented such that the +Y^ plane is. taken 
as the model plane of symmetry (Ref. MSEC 14A00305). The {+) is con- 
sidered 0 deg, the keel side, and (-)Yg as the 180 deg reference, the lee side 
for the pressure loading. In this orientation the skirt holddown posts are 
located at 60 and 120 deg, see Fig. 3-2. Another SRB model has been con- 
structed in which the plane of symmetry has been rotated 90 deg so that the 
+Zg plane is taken as the plane of symmetry. This model is designated as 
the "rotated model." In this configuration the holddown posts are located at 
30 and 150 deg. The 0 deg reference in this model is also considered the 
keel side. This configuration is shown in Fig. 3-3.^ The total number of 
rows and columns has been held constant for both models with the only 
changes being made in the nodal column locations. 


3-3 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



rH 


LMSC-HREC TR D497054 



Fig.3-2a - Original Model Branches 1 and 2: Grid Layout 
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Note; Heavy lines denote stiffeners. 


Fig. 3 -2b - Original Model Branches 3 and 4: Grid Layout 
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Fig. 3 -3a - Rotated Model Branches 1 and 2: Grid Layout 
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Fig.3-3b - Rotated Model, Branches 3 and 4: Grid Layout 
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The reference surface for each branch in the model was assumed to 
be the outside skin surface. This was chosen to simplify the data input in 
the skirt branch and reduce the possibility of having input data errors. The 
outside skin surface of the conical skirt is a constant radius at any given 
longitudinal station; while the inside surface is stepped for the changing thick- 
ness in both the longitudinal and circumferential directions. 

/ 

Tables 3-1 and 3-2 are listings of the WALL subroutine for the original 
and rotated models, respectively. The variation of wall thickness versus X 
and Y coordinates for Branches 2 and 4 and the smeared stiffener data for 
Branch 4 are given. Tables 3-3 and 3-4 are the corresponding additional 
STAGS C input data listings for a linear analysis for both models, respec- 
tively. Refer to the STAGS C user's manual for interpretation of the data 
listings. Note the ordered sequence of the branches and the interconnection 
of the node numbering sequences from one branch to another. This is the 
most efficient sequence for model accuracy and computer run time. The re- 
sults of the parameter and grid spacing studies reported in Section 2 were 
used in determining the nodal locations for the model. 

3.3 MODEL RUN TIMES 

Typical run times on the NASA-MSFC Univac 1108 system for a linear 
analysis using this model was 20 minutes CPU and 13 minutes I/O for a total 
time charge of 33 minutes. A bifurcation buckling analysis to determine the 
first eigenvalue and mode shape took 38 minutes CPU and 40 minutes I/O. 

3.4 MATERIAL PROPERTIES 

The material properties used in the model and the allowable stress 
values for checking the results are given in Table 3-5 for both the D6AC 
steel motor case material (Branches 1, 2 and 3) and the 2219 aluminum skirt 
structure (Branch 4). 
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Table 3-1 

WALL SUBROUTINE LISTING, ORIGINAL MODEL (+Y_) 

— & 


SFOR.JHS HALL, /ALL 
HSA -E 3 - 10 /I&/ 7 &~atU! 

-i*7 

C - SET JF> SECTIONAL-PROPEFLT-KS-tW-BA*-A--STWW-£HtNM- 
C C 3 L = COLUHN LOCATIONS AL 0 N 5 Cl RCUHFERENCE 

C ROW -i-ROW-l 

C THICK = THICKNESS EOS SECTIONS ACCORDING TO <ROW,COLI 

C ISIF = NO. OF -STI-FCNERS-IN— COLUHN -SECTION 

C K S I F = S&SINNING AMD CHOIRS RONS FOR SHEARED STIFFENER’S 


r - 
C 


DIMENSION I$TF(8l,XSTFm 

■ - 01 HE N SI ON COLtl7>, »0«M4J 

OUENSIQN THICK(16,14) 

- - DATA COL/O.D, 7.5, 2 2. 5, 42.0,54.0,60.-0,^8. 0,7».-3,90.0-, 102 . 0 , 11 4-^0i 


i. 

— -DATA 


120. 0,126. 0,138. 0,157. 5,172. 5, 180. 3/ 
R0W/3.a,7.72J»B.2-5,18-.-3O5, 


1 , SB. 0,83. 39, 83. 9, 90. 9S6/ 

- OATA THICK/16 + 1.375, 2*0. 9, 2*0. 9475,-2*l-^I75r4-*O.WS- r 
2, 0.9875, 2*0. 90, 


3 2*3.4 25*2*0.6,2*1.375,4*a.6,2+1.375,.6tJTrO.S»2*3.4-2Sr—— 
42*D.4S25f2+3.55»2*1.375»4+G.55, 2*1. 375, .55, .55,2*0.4625, 


5>i*3.5, 2*1.375, <* + 0«5 t 2*l»375j4*»53» — — 

5 2* 3. 4525, 0.4695, 0.5, 2*1. 375, D. 5, 2*0. 4625, .5, 2+1. 375,. 50, .4 68 5, 2* 

72 + 3.4 2 5,3. 437, C.S, 2 + 1.37 5,0.5, 2, *0. 425,3’. 5, 2*1. 375..S ,.4 37 ,2*. 4 25, 

42*3. 4625, -3 .4695, 0.5, 2*1.375,0. 5 ,2*0-r4-525,-3-.5 ,2 *l-»-J75 ,0.-5, 3.4-68-5, — 
42+3.4525, 

- 94*3.5, 2 + 1.375, 4 + Q * -5, 2 ♦ 1 » 3 7 5 , 4 * 0 . 5 , ■ “ " 

13*3. 4525, 3. 5, 2+1. 375, 0.5, 2*3. 4625, 0.5, 2*1.375,0. 5, 3 *3. 4525, 

1 3*3.425,0.5,2*1.375,0.5,2*0. 4 25 ,0.5,2*! .37 5 »0.5,-3-*D.425, 

23+3.4525,0. 5, 2+1.375, 3. 5, 2*3. 4625, 0.5, 2*1.375, 3. 5, 3*0. 4625, 

J-i + 3. 5, 2*1. 375, 4*0. 5, 2*1.375, 4*Q. 5-, 

94*3.5,2*1.375,4*0.5,2*1.375,4*0.5/ 

3ATA-ISIF/-1, 3,3, 2,0, 0,2,2-/ 

DATA XSTF/7.723, 18. 006,24. 006,51.481, 57.481*83.390 ./ 

DIMENSION DRQwU 1 > ,OTMXCM 1 l > * -- — -•> — 

OATA DRDH/93.3, 91.086,92. 174, 101.872, 113.478, 125.185, 

1 128.0,133.0,131.0, 132. 0,133. 171/ 

DATA 3 THICK/6 *0.362, 0.42,0. 6,0. 9., 1.75,3.0/ 


C CHANGE THIS CARD FOR EXECUTION OF STAGS FOR 4 BRANCHES 

1 F < I 3RNCH . £3 • 4 ) GO TO 75 

- IF (I3RNCH.EQ.2I GO - T 0 76- 

RETU-'N 

74-. CONTINUE 

I C F 9 -3 / 

30 7 1 1 = 1 ,11 — f 

TLin.'O THICK (XI 

r - t ! X-OTOi)<IU .LT.O.OOU-SO-TO— — 

71 CONTINUE 

--7 2 2 = -TLU 1/2. -3 : 

•K4(lJ=j.3EB 

iVl(U=3.iL5 “4-= 




3-9 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 

i 1 



r> r» C 5 


LMSC-HREC TR D497054 


Table 3-1 (Continued) 


U?U1»=0.33 

53U)=0.i1E8 

rEi3tn=o.o 

SH 33 tn =»-0 1 

LMS=1 

- LSI } S -0 

ISTFF=0 

I»U=0 

R4DA=3.3 

- - • return 

15 CD Ml I M UE 

10-3=3 — 

l»T5-l 

- ■. LSrRS = 0 " — — 1 — 

.I D ?w:0 

^IlU^ftl) COL 

rf?ITE(6,621 ROW 
52 F DR^AT l *- ROM * » 16F7-.-4I 

rfRlT£«S,£>OI < < THICK! J.IJ. J = 

50 FORMAT! ‘ THICK *Tlt>FT-*T 

5 CONTINUE 

'JR0K=1>* 


NC0L=17 

00 1 i=i,mcol-- : 

I F ( tr.SE.COLtm-rAHO-sfY-itfTreOtrf-t+l H J 30 TO 3- 

I CONTINUE 

3 CONTINUE " 


IC3L IS THE COLUMN NO. OF V 

DO 2 I=1,NR0H 

I =i T * -t 

IFUK.BE. ROW «-I-Hv«*9v* frrtTTriHHtf M-H+6 ~ 

C 0 *JT I NU£ 

• I R 0 V I S-Frt E-RON-SE CTION NOi- BF- K- ! NOT - NC CCS5 ** 
EXl=13.5Eb 
-£X1 = 10.5£5 ■ 

XN1 = D. 5 
El =EK1 
Di =XN1 


■ s=- E X i / HHHrH — 

2 £ T 2 =0 . Q 

IC 1 = 1 COL 

IC2=IC0L-l ’ „ 

IF (IC1.EQ.1) IC-2=I 

IFIIC1.E3.17)1C1=16 

- AT ;0.5 + !THICK-H-e-H4fl0yl t Till CK 1 1 Cl 1 1 ROM H 

Ar(COL(ICOLU>-COL!ICOLn/!COL!ICOL*l)-COl^«L-in 

- 3 S( C 0 L!IC 3 L»-C 0 LUC 0 t-HW« 0 t*mfcM 4 - 6 Ofc«-Mfc-W-- 

AT = 1 . 0 *!A»THICKUCl ? iaO«l*a»THIC l!iIC 2 ,iaOUM 

ThI^CONS'LETES CALC. FOR HONOCOQUE SHELL 
-NO* ADO IN-SH5-ftR& 0 STIFFE NC -W “ 

HOJE THIS CARO UP 


original pact? to 
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Table 3-1 (Concluded) 


c 

iri=i 

J<1=0 

m =1 - m 3 

— — 511=0.-52424 

51 = 2 m S933E6 

— E2l =-l .730 

E?21=-2.453 


01=3.3— 

rilli:ll 

. EX3<1I=EK1 

EY3(11=EU 

Jnui:<Nl ^ 

63(1 1-6 

ZEI3<11=0.0 

^H33CU=0.D 

- - - - IF ( ( X.GE.XSTF1122-.414O****l-2*XSTF-12>42g0 TO 10 
lF<(X.G£.XSTF13)).A4D.iX.LE.XSTFm))S0 TO 13 

IF ( (X.GE.XSTF 1 5) 7TrAttO»4X-»fcE-» X5TFt 61 - l -l B O - TO liT 

I S I FF- 3 

’SO TO -5 0 

13 JC=ICOL 

IF ( IC0L.5T.9) IC=17— TOOL 

C IC LOCATES COL UHN— SEGT-I-O-M— MLIMBES 

I0=1CDL-1 

IF (TO. 37.9)10 = 1-7— 10 


I p I ( 2C3L.E0. 11.07. ( ICOL.EO.17) ) 30 TO 20 

-- - i'( ( I5TF (ICJ »ISTFHD> V.tO.Q ) — SO -TO -2-1 2 

01 =1 .0*1 COL (ICO L» 11 -COLtZCOl-in /(ISTMIC)*ISTFIIDI 1 

• 50 TO 21 - 

20 01= (C0L(2)-C0LC11 1/1.0 *ISTFC11 

- — 21- COITIlilE T 

1511=1 

- • ■ ' - ISH1(1>=2B - — ■*» " i. — * *J- . dm rn ' d dSddl d dlil- d d — du 

isiacn=i 


SMS B 1 1 1 = 01 - — 

51416(11=3.3 

I F i 43S<r- 34.001 .LT. 0.051 01= Q.S*f 93.30- 7B. 001/1. 00 

r F ( A B S < ¥ 0.0D 1-.L-T. 0.054— 01= — 0-.5*-l — 95*-OD—- a4»0&F/4.0Q- 

I F ( ASS<¥- 96. 031.LT. 3. 351 01= 0.5*1102.03- 93.301/1.00 
-- - I F ( A3 S 1 ¥ -1 02*03) .L-T •.0.05 1- Dl=— 0.5*14-l-4-.00-— 96-*3B4-/-l-o5-Q— 
I F ( A351V-138.33) .LT .3.05) 01= 0 . 5*( l 4 7.75-1 2 5 . 33 1 / 1 .75 
JF( . A B S < ¥ - 14 7 .7 S > *UL^Ua5)-01— 0.5*4-1.53*25^433 .3 3X41 *5-0— 


I F C ASSCr-157.531 .LT.3.05) 01=0.5*1165.33-147.751/1 .25 
-- ~ I F ( A 3 S ( Y -165.001 .LT-.3.05) Dl=- 0.5* ( 1 72 .53-1 57. 50 1 / 1 *00 - 

I F ( ADS(¥-172. 50) .LT.3.05) 01= 0 .5 * (1 8 D . 33-1 6 5 . 00 1 / 1 .00 

53 COiriNOi ~ 

C 50 W4ITE(£>,S3) it, r, 101,102, AT, IC, 10, 01 

IF ( J1.LT.0.01>-ISTFF=3 — 


53 FOIIAJI* *,2(E14.5,2X),2I6,2X,E14.6,2I5,2X,F14.5) 
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Table 3-2 

WALL SUBROUTINE LISTING, ROTATED MODEL (+Z B ) 


UFOK.UtfS WALL i MALL 

USA C3 -10/16/76-(r3;T?1-2l— ftr.T) 

-47 . . — 

C SET UP SECT I ONAL” PK0PERT1E 5~TS " 0'XT7C~'5 Y’A’TE HE N TS 

C COL » COLUMN LOCATIONS ALONG CIRCUMFERENCE 

-C R 0 T ROTi SrCTTUN~LTJC AT IONS KCO snJT.1t N OT m" (Nof HITE'S S' A R I L Y ROW NO*} 

C THICK * THICKNESS FOR SECTIONS ACCORDING T« (RO*,COH 

C " ’ ISTF - NO. OF STIFFENERS " JN 'COLOUR 'SECT fOR — 

C XSfF » BUG I NN 1 NG and ending ROWS FOR SMEARED STIFFENERS 

c 

DIMENSION ISTFUOJ (X5TF(6) / 

■ DIMENSION C0LI2I 1 ,K0MI4) 

DIMENSION th i cm 20 >T 4 1 

data COL/Q.O, 12. 0^24.0,' '30.0, 3ft .'(U48 ;0 ,67~*50 »'b2 0 ,97 

II 12, 5, 132. 0, |38. 0, HR. 0,150. 0,156.0, 162. 0,168.0,174, 0.1 80.0/ 

“DATA ROwVOiC, 7. 723', B‘, 25, ffi.u'QBVrBiTB , 24.006 ,24 *5 ,Rl ^a'l ,'52,o7£7»481 
I ,58, 0,63.34,83. 9, 90. 954/ 

" — DATA THICK/ 20*1.376,' ' 

22 * .98/5, 2 . 1 .37S»2«.9a75,5*,9*3*.9875.2*1.37s,‘<*.9875» 

32*. 6, 2*1. 3/S, 2*. 6, S*. <(25, "3»V<>J2*T737,S';‘»*T{T ' 

4 7 *. 55, 2*1. 375, 2*. 55, 5*. 46 2b, 3».S5, 2*1. 375,4*. 55, 

S t * « 5 ■ 2 » i « 3 7 5”, 1 0 * ‘. ST 2 « i 73 7 5" . 9 * . 57 

6, 46 25. » 6. 2* I. 3 75,, 5, 7 ..4o25, 2* .5, 2* I. 375, 2*. 5. 2*. 4625, 

- - 7 . 4 25 , . 5 , 2 * 1 . 3 7 <5 » . S . . 4 3 7 ,'5 * . 42 5 174377 2* 757 2 *T73 7572*7*72 *T4'2 57 

D. 4625, .5 ,2* I .375,. 5, , 46b 5 , S* , 4 625 , . 4 685 , 2* • 5 , 2 • 1 . 37 5 , 2* • 5 , 2* » 4625 , 
92««S.2*l»37S, 10*. 5.2*(, 375,4*. 57 ‘ 

1. 4025. . 5.2*1.375. .5 ,7*. 4625,2*. 5, 2* 1.375, 2*.5,2*.46ZS, 

1 . 425 , .5,2* 1 . 375, .577* .’4"2s7 2*7 5 .2*17375', 2*.5", 2* 7425 " 

2. 46 25.. 5 .2*1 .375, .5,7*. 4 o25, 2*. 5, 2* l .375, 2*. 5.2*. 9625, 

" ~ ‘ -52»jS,_M .375, |£J*. 5.2*1 .376.4*. 5, 

4/*. 5,2*1*375, 10*. 5. 2*1. 375,9*. 5 / 

DATA I5TF/2*Z,2*0,2,J,2,4*I,3j2*l 2*079 *1 / 

DATA X5TF/7. 723.lU.Ou6, 24. 0U6.51 ,481 ,57*981 ,83*390 / 

DIMENSION OROftdl f.DTnlCKt in ~ 

OAt-A OROrt/Vo. 0,91. 0B6, 9?; 179. 101 .872,113*978,125*185, 

1 12&.0, 130. 0,131.0, 132. 0,133, 171/ 

DATA DrrtJCK/6*0. 362,0. 42,0«6, G»4, 1.7S, 3.0/ - < 

C CHANGE THIS CARO FOR EXECUTION OF STA GS FOR 9 BRANCHES 

JF< IBHNCH.EQ.4l 60 TO 75 

- ■ IFt I URNCH.Eg.2l GO TO 76 

RETURN 

’ "76 continue L 

ICFu»3 

00 71 r-Till 

TLU l*DThtCKU )'■ 

IM (X-OROv.[ H r f LT,"OiOOU _ GO T0'T2 —— — — — 

71 CONTINUE ’ 

* ''72 ZWLUI/2,0 : 

Ex3(lJ*U.3£0 __ 

EY3 ( l I *0* j£6' — 

U2 ! ( 1) “0. 33 

gj( ifo.utu — ■" ' ; 
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Table 3-2 (Continued) 


££T3 { l l “0*0 

Rnoatt )*o»o 

LAtS=t 

— - lst rs*o 

tSTFF-O 

IPKffiO 

khoa-o.q 

7b CUNTINUE 

|CFB=3 ■' 

LATS-l 

L'STKS “0 — — 

IPNrt«Q 

u t u.GT*ff¥Di",0Ri"nvsr^iBr)''.^t)-txr5" — 

WKITE.16,61) COL 

‘61 FORMAT t * E0(_* ,l|T7«'l ) 

*Hl Th (6,6?) ROW 

■'T2‘>oriMAn rftDiifjTiFrmi ______ 

WRITE (6 ,60) ( ( THICK ( 0,11 ,J« 1,20 > ,1*>I , ) ** ) 

60 FORMAT! 1 THlcK r i 20 r 6 . 3 ‘J 

S COMT 1 RUE 

HKOfl.H-- ‘ 

M 1 0 L “ 2 ) 

do i _ r*i»NcoL 

ICOLbI 

if( ( y . uE , col ( t ) ) .ANortT.LT^coirrryi ' n ~t 6 t)~nr~3 

1 CONTINUE 

'“3 CONTINUE " — : 1 — 

c • icol is The column no. of y 

00 ‘ 2 iPljNROW _____ 

I K 0 A • 1 

2 1 F ( ( X . GE . HOW ) D ) i'A N CT (TCTCTTR 0 WTr+TTlT<rCTTD '4~ 

•i CONTINUt 

C IROw IS the row SECT I 0N~ND • 0 F “X ~TTT0T~ HECgS SA H TU ' Y ' T M E~R CTW~WKFE R V 

EX l n 1 U. SE6 

. LY 1 « 1 U ■. 5 fc 6' ‘ 1 

A N i « U • 3 

1 U1“XN) 

G* tXl/(2.*(l,**Nl»| ’ — - 

2E I 2»U*l) 

IC1 „, C0L - 

IC2*JC0L-| 

■ — iruci,e«.nic2 i, i 

lf('Cl,t«,2l| I C 1 *20 

C Al*=O.S»(ThlCK(ICI ,lROW)FTHtCK11cr,TR0Wn 

A=(COL(lCOL*I J -COL (i COL))/ (COL ( I CQL+ i ) -COL ( I COL* I ) ) 

— b*(C 0 L{ IC 0 L)-COLriC 0 C--tT) 7 TCOr(TC 0 L*I l-COVTTCOrSTT) 

AT-l.U»<AeTHICK(lCl , I ROW ) +B »TH I CK ( I C2 , J RO«k ) ) 

X -“AT /2 , •- - 

This COMPLETES CALC, FOR MOnOCOQUE SHELL 

r.CiA J UU Ifi SMEAKtp STIFFENERS 

,:uyc this card up 


c 

c 


3-13 


LOCKHEED -HUNTSVILLE RESEARCH, & ENGINEERING CENTER 



LMSC-HREC TR D497054 


Table 3-2 (Concluded) 


' JKl»0 ” “ ' 
Al'UHO 

SI I « 0 •B 2 HSH 

SI-2.6903E6 

TEZli-517W 

E2/l=-2.4U0 


01-0. O' ' . 

fL(l)«AT 

CX 3 t -j ) «EX1 ' 

U2I(11«XN| 

gj<i >«g 

■“ ZET3U1-CUQ 

Krt03< 1 )» 0.0 

■ Tf'UjUGLVxS'TFd J ) .AND. (X,L£.XS'tF(2) ) )GO TO 10 
IF < U.GE.XSTFO) I .AND.lX.LE.XSTFlR) ) >50 TO 10 
IF t (X«G£»X5TF (5) > * AND #"( X »L"£ • XSTF’i' 6 I ) ) GO TO IQ 
I STFF a C 

to ro so ' ■ ' “ 

10 IC-ICOL 

IbTFF‘1 

ic locates column section number 

IL>»IC0L-I • 

IMUCOl.E(I# 1).OR,UCOL.EQ.21)) GO TO 20 
I F 1 1 ISTFllO'HSTFl 1DJ ) ,EQ,0| ~50'T0~ZT 

Dl«l .DMCOLllCOL-M >-C0L< ICOL-l ) > / 1 1 STF 1 1 C UlSTF U 0 )) 

00 TO 21 - 

20 jr UCOL.tG.J) 0l“(C0L12)-C0L( 1 ) J/1STFC1 ) 

1 FUC 0 l.EU. 2 J >Dr*(C0L(2n-C0L<2ffTJ7 _ TSTFrinn 

21 CONTINUE 

Nsri-t — 

! S MN l I ) »2 4 
1SM0U >-l" ” 

SMSP ( 1 J -D 1 

■ •“ ■ 5 tl AM G I 1 > a 0.0 

SO CONTINUE 

' IVRITEU.63) X't T * I'C I » I C27AT i IC’iTD 7D 1 * 

1FCD1 .LTfO.Ol) ISTFF-O 

1\J - r 0 K H ATT * -» *2 ( E r'nrs TZ X TTZTi . 2 X , E l STS! 2 iv."2 XiEiH.il 

RETURN 
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Table 3-3 

STAGS C INPUT DATA LISTING, ORIGINAL MODEL (+Y B ) 


$ T A3S OITA. 


* * * * * * - 

5 T t -&s — £— e EW-tt-s a a *oo e-l— • 

0 4 *i 3 Q Q 

. 3 21- — u — 2 l- 7 21 38 21 

1 b 3 D 3 *' 3 8 

2 1- 0 —3 1 3 8 

3 10 3 13 1 

H 4 37 7-2 24— —O * 

13 2 1 

3 3 4 1 

5 I D 3 

r 3 -214 . 35 .4- WO r Wo — 


44-72 44- 2 2- 

84.24 93.24 1D2.25 

4-2 S. Jb 1-3-1 o-27- 14 2. 27-- 

169. 2B 174.29 194.29 

23 9.8 2 212-S-9 — 



132. 133.171 

7,5 -2 2.5 - - 92o0 
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Table 3-3 (Continued) 


STAGS INPUT DATA. 


PAGE 2 


* * 

* 

* 

+ 

4 

* 

* ' 

, a i xa-n 

6**0 9 3 • 

147.75 ' 157 

0 Mrr 

.5 115 

tj rircrg ht 

.□ 172.5 180 

.0 




0 1 







i 

• — 1 5 

2.523 -3. 

“ * 3 3U.~5ro 

139 3,2996 -2. 

1 r.tt - - 

5392 —2 . 

523 0.139 -0. 

2996 2. 

5092 

1 2.-577 13i 

1 11 

■mrx n» 

1 21 

TJl CT3 






0 — -- T- 

. 1 3 

—l~~ 1 







- ' - 5 — 1 

.0 

’O' * 3 — 

13.49 

.0 

130. 

73. 




— — *t 
0.0 

1.311 

2.622 

5.245 

7.867 

9.178 

nn- 

10.49 
-66 - 

7-8^- 

■ 0.0 ' 
8 4. 

7 *'5 

90. 

- — Zc.tr - 

96. 

■ r He. U* * — 

102.* 

114. 

120. 

126. 

138. 

147.75 

6 4 

t57'.5" 
6 4 

1 6*5 . 

ItZ . a — ■ 





o i— 

1 0 

4 3 

30. E6 

11. £6 

.28 

0 

- * 


5 *625’ 

,8.571 

1-» 42 5- 

57.56 

3.333 

.0 ■ 

1.94 - 

.0 

2.095 


0 1‘ 

— i 2i 

1 1 







— * t — 

3 a 








.52 

—.25 
33. E6 

.33 . 

30. E6 

11. E6 

r 

- 


.0 

□ — *o 

90.956 

.0 

180. 


104.1 



— t- --I 

.0 

1.5 

3.9 

7.25 

7.723 

P 1 c 

8.25 

XI % 

12.36 

17.5 

2-4 .5 

1 a a G06 

30.87 

O'. 3- CUI3 

37.74 44.61 

SO. 9 

51.481 

a.ot 

52. 

70*4-? 

5 4. 

54.431 

S3. 39 

83.9 66. 

35.536 

87.0 

C_u 

90.956 
60^ — 

Sfr , 

7-S-.-- 

• 3 

84. 

93. 

95. 

102. 


120. 

126. 

138. 

147-. 75 

6 4 

1 5-r-*"5 

3 4 

' ibr**- 






-6874. 

i i 

12 16 






j. g 

-7.389 

5.297 

•’ 3 . 3 4 9 

.087 

-7.969 

2.692 

■ __n. 

-.639 
-H »357 

.21 6 

— 6*-9 2 4 

1 3 

— Cro-* rc c 

4 3 

Tq.Es 

4 . E 6 




.,314 

— - — - —7* 5 8 5 

6.726 

5.465 — 

84.243 

19.469 

25.082 

.275 

-1.37 

-4.362 


-0.86? 

6.392 

-13.53? .439 

—9 « S 2 ? 

3.253 -.539 

.216 

1 1 * 2 7fr 

—I *fr2 .895 
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Table 3-3 (Continued) 


STAGS INPUT DATA. P* S£ $ 


i 3— 

13.963 

— s a — 

285.993 

1,0.66 — 

7% *5 OD 

* n re 



59.830 

t\ , ££ 

226. 

-1.474 

4.363 


-2.73 

1.25 

-2 .f 3 

-1.25 

-2.70 

, 1 ?*l 

0.0 

.0 

-.5 

— -2 . 299 

'0.3 

J"1 O 79 — ' 

1 3 

— ■ — r99*J-" 

4 3 

ID. 65 

4.E6 

-2*7^ 

_3 

=4-. 

UUO- 

-f -J— — 

2.713 
V 3 

1.172 

4 3 — 

.994 

10.ES— 

.0 

■4-.E6- 

.421 

.0 

-2.174 




25.42 


85.511 


-i - -a- — 

16.715 

-0 3 - 

62.267 

10.66 

9.735 

— 4.66 — 

.0 

34.3 

.0 

3.343 

1 7.717 

Q J.- — 

74.155 

1 0 « £ 6 

9.229 

■ 1 *? rt O"’ 

.0 

36.9 

.0 

- 3.543 

• 4.973 

t n 

— □ -3 

26.185 

. 2.5B9 
i.n _F£ 

.0 

4 .E 6 

10.3 

.0 

3.976 
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Table 3-3 (Continued) 


STUBS INPUT DATA* PA BE % 
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Table 3-3 (Continued) 


STA6S INPUT DATA, PA5E 5 
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Table 3-4 

STAGS C INPUT DATA LISTING, ROTATED MODEL (+Z B ) 


WELT • 0 I L TPFS.ARTSKT 
fccTQ07-KLTB90— 10V-l6w«*ll2M-5--t“.e) 


STAGS C SRB ROTATED TO DEG 


- 0 --•*»“ 

12 21 

n 2 i 

0 V 

7 21 

36 2 1 





1 6 """ 

2 l 

3 1 

4 6 

0 u 

0 0 

3 0 

1 0 

1 

1 





0 "" "0 
HO 72 

1 ' 0 ' 
25 0 

1 





1 3 

1 3 

2 l 

3 1 







3 3' 

1 . 

- 5 - i -- 

.0 

*T “ V 

0 , 

U 






0 0 
215.35 

.0 

iao. 

73, 




- 1 * I 

Q.C 

2 « 5 

5«0 

16. 

26.22 

39.22 

11.72 

11.22 

16. 72 
111.25 

49 ,“ZZ““ 
US, 75 

~51T22 

ua.oo 

~5T.2‘3 

120.26 

~7'1'."Z7 

123,01 

BITJfT" 

125.76 

rsrpr 

131.27 

102.25 

112.27 

1S3'. 2 ? ~ 
191,2V 

158 . 77 " 
199,30 

161,52 

201.80 

"m.' 2 e 

201.30 

'I66.7B— 

205,68 

16 9, ZB 
207.06 

T71TZ9' 

209.82 

. “”T84V2?“ 
212.58 

2 1 3 » y 6 

« 0 

215.35“ 

12 . 

21 . 

30, 

36, 

18 . 

67.5 

82»5 

VO* 
l B 6 • 

97 »5 
162 . 

105. 

168. 

11215' 

171. 

132. 

160, 

136. 

Til. 

150, 

2 H 

a i 

'6 ‘ 4 







1 0 
• 5H6 

1 " 0 “ 
2. ICS 

30 • fc .6 

11*66 

2,105 

* 2 & 

- 1 * 0^6 

-2.105 

.516 

-2.105 

6*7 

i o 

— 1 *~R 1 5 
4 Q 

V.BV9 

3G.E6 

•0 
1 1 ,E 6 

SET519 

.28 



-1.825 

3 • 65 
i*VV3 

* n -* 

1 o S 2 5 

5*6V2 

-*■ i n 

" »r 
1.621 

1 • 8 2. 5 
•0 . 

. r~ 

.0137 

~T»B25 
• 0 

3*05~ 

1*975 


1,481 <0 


“ «»/l .16*** fS737 iO ~ — TOTS f 0 ♦TRB' 

3 l l 21 


- — 3“ 

1 

8 —■ 

11 - 
20 “ 
28 

* 1 - 21 — 
1 21 





- 

2 

2 

I 2\ 

\ 21 



* 



1 

3 6 

i 2 1 






0 

i 

1 1 






'“”3 ~ 

0 







i 

i 

- .26 







*B 2 

30 *E 6 

V33 30.E6 

1 1 * £6 

• 0 

. 2 * 


7 

t 

0 Q 






TU*U 

T33~ «" 1 7 t 

'■'O'.B EBTOro 

73.0 




-1 

-i 







7U»U 

9 1 ,tJS 6 

92“oT71 10T.S7Z 

TrniTa 

126 • 

1 26 9 

“"130.“ 

131 . 

137. 

133.17J 






« u 

* l 2 * • 

“21. ”30, 

36 * 

' 10 . 

67 • & 

■ Bars - 


VO. 

97*S 

105. 112.5 

132. 

138. 

111 . 

150. 

* • 

1 6 x * 

“TifS, 17T, 

1 UUl 




6 

1 

3 1 






' 0 

l ” 

- — - 



... 

— 

— 

— 
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Table 3-4 (Continued) 


1 0 

4 0 • 

30. E6 

11.06 

• 28 




•' 2 .525 

12*577 

' -i'i'3? 
13.5321 

' .2996 
13.237 

-2,5(392 

2,3068 

-2*5X1 

21,825 

,1'39” 

-.2996 2.5092 

- l - ri - 

0 1 

1 21” 
1 1 







— j jt 

5 1 

0 0 







-1 -l 

TO. 4 9 

.0 

180. 

J 

73, 




0*0 

.0 

i *'3i r " 
12. 

2.622 

24, 

Tf» 2 4 5 
30. 

7*tt67 
36 • 

9,t7e 

48. 

10.49 

67.6 

82.5 

■ ""9’u» 

156. 

n ;s— 

162. 

105, 

160, 

112.5 

174, 

13 Ti 
180, 

138. 

144, 

150, 

6 ~ * "4" 

0 t 

4 ~ 







j-. 0 ---- 

5*625 

' “4 0" 

1 .425 

£0 • E6 
5,62b 

U".E6 

-1.425 

,28 

5,625 

• 0 

.0 

.0 

S'* Si? 1 
l 7 

57.56 
1 21 

3,033 

*0 

1.94 

.0 

2,095 


- ‘O' ■ i ' 

3 0 

"1" 1“ 







~ r - i 
a 52 

- " -.26" 
3Q.E6 

,33 

30. E6 

11. E6 




6 T 

.0 

"O' D~ 

90.956 

,0 

180, 

73, 

104.1 



■-T“ -l ’ 

• a 

1 • ^ 

3.9 

7.25 

7.723 

8,25 

12,86 

17.5 

" 10.006. 
30.87 

13.5 

37.74 

20. b 21.006 

44.61 50.9 

2 1 9 5 
51 *481 

23.5 

52. 

2«7qq6 

54* 

24 «S 
54.481 

"" 55. 
83.39 

”57. 

83*9 

5774 BT 

86, 

58. 

86.536 

63.95 

87.0 

70.43 

90,956 

7679 l 

82.9 

■ .0 

90. 

1 2, ” 

97*5 

247 

105, 

30. 

112.5 

36. 

132. 

.48 , 
138. 

“ *67,5 

144, 

82.5 

150. 

•• 156. "■ 
6 4 

"162. " 
3 4 

’ "168, 

174. " 

lUO. 




“ c • ' I ~ 

1 0 

4 0 

10. E6 

4.E6 





‘ -7*089"" 

6*924 

"'5.297'" 

68.722 

" -8". 8 49" 

17.3 

— '.06 7 ' 
19,61 

-7*969 
• 326 

'. 2 .692' 
“1,37 

- * 6 39 
-H .057 

• 214 

' i - • o — 

-7 *586 

4 li" 

5.4*5 

10.E6 

-9.346 

"4 > E 6 ” 

.255 

*0.466 

2.860 

- « 639 

*216 

— So' 726 

1 0 

04.740— 
4 0 

T0T469" 

10.E6 

Z^*'0 8 2" 

H • £6 

*27 b 

-1,37 

~-*r. 06'2 


- -8*662 
1 W76 

6.092 

142.896 

— 5S"t 0~* 5 a 2“ 
33.596 

.'408' 

41.681 

*9 • 6'22 
1.477 

3. 2 50 
-27393 

"■ ♦ 63 V 

• 2 1 6 

- -J - - o - 

13*960 

o -- 0 
265.970 

'" " 1D.E6 
74.5QQ 

4 1 E 6 
68 • 83q 

226. 
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Table 3-4 (Continued) 
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Table 3-4 (Concluded) 
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Table 3-5 

MATERIAL PROPERTIES 


Property 

Material 

D6AC Steel 

2219-T87 Aluminum 

E (I0 b ), (psi) 

30 

10 

p, (lb/in 3 ) 

0.28 

0.102 

F tu , (ksi) 

195 

63 

F t y’ < ksi ) 

180 ' 

51 


3-26 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 










LMSC-HREC TR D497054 


4. SRB CAVITY COLLAPSE LOAD ANALYSIS 

4.1 SRB CAVITY COLLAPSE LOADS 

The SRB mathematical model was analyzed for the nominal conditions 
of 85 ft/sec vertical velocity at water impact as specified in Revision D of 
the SRB loads document (Ref. 1). The distribution of the external pressure 
on the components of the SRB vehicle during cavity collapse are shown in 
Figs. 4- la and 4- lb. These figures are reprints of data presented in Ref. 1 
and show the pressure as a function of the wetted length along the keel side 
(< f> = 0) and lee side ($ = 180) of the vehicle. The circumferential variations 
of these pressures on each component are described in functional form also. 
The keel and lee meridians may fall along any surface meridian of the SRB 
vehicle. For these analyses the lee and keel meridians were assumed to 
fall on cross sectional planes of symmetry (FYg or Internal skirt 

and bulkhead pressure loads were determined from the hydrostatic pres- 
sures, P^ and Pj given in Ref. 1, and the functional circumferential distri- 
butions as specified in Fig. 4-1. The external circumferential cavity collapse 
distribution is shown in Fig. 4-lb. A constant internal motor case pressure 
of p = -1.65 psig, as determined from Ref. 1, was used in the analysis. This 
negative internal pressure is additive to the effects of the external pressure. 

The SRB cavity collapse load for the 85/32/5 (VV/VH/9) condition is 
being analyzed. This case was determined to be the maximum pressure 
differential for the nominal 85 ft/sec vertical impact velocity condition. 

This determination was made by comparing the external pressure to the 
internal skirt hydrostatic pressure for each condition. Since all of the load 
cases have the same distribution and vary in magnitude of pressure levels 
only, this rough test for maximum load case is valid. The pressure differ- 
ential for this method is Ap = 148.5 for the 85/32/5 case. The load values 
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Cavity Collapse 

SRM Case and SRB Aft Skirt Pressures 



INTERNAL PRESSURES 


SRM CASE 


SRB AFT SKIRT 


P - ^ULLAGE' UNIFORM 


P » Ph + ^(1 * cos $)(Pi - Ph) 
0 < X < Lp, -130° < 0 £ 180° 


external pressures 

CIRCUMFERENTIAL DISTRIBUTION 

P = P KEEL + K C (P LEE ’ P KEEL> * 0 <X < Ly^, -180° < 0 < 180° 
p =-P KEEL» Lw,L-X - E W) K t - 0VJ - 0 < 0W 


-1 f 2C X - Lw , L) 

V SEE FIGURE V.C.III.2; «„ = cos ' Tw) 

LONGITUDINAL DISTRIBUTION • 

SEE FIGURES FOR PkEEL & P LEE vs * STATION 


L 


Fig. 4 -la - Cavity Collapse Pressure Distributions — SRM Case and Aft 
Skirt Pressures (Reprint from Ref. 1) 
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BULKHEAD PRESSURE 



p ~ ^ULLAGE UNIFORM INTERNALLY 


EXTERNAL 

p ='?i + %<Ph - Pi) (I - cos 0),, -180° £ (S i 180° 
UNIFORM LONGITUDINALLY 
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00 
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External Circumferential Distribution 


Fig. 4- lb - Cavity Collapse Pressure Distributions — SRM Aft Bulkhead 
(Reprint from Ref. 1) 
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are taken from Ref. 1. The pressure factor, ratio of the maximum pressure to 
the maximum value given, for the nominal distribution as shown in Ref. 1, is 
1.176 for the 85/32/5 case. This factor is to be applied to the external pres- 
sure distribution only. Note no design factors have been applied to the loads 
in the analyses. 

The nominal pressure distribution is to be modified by shifting the maxi- 
mum pressure peak location for the lee side over a range of +l/4 SRM diameter. 
This shift will account for uncertainties in the prediction of vehicle attitude and 
ocean conditions at water impact and cavity collapse. The fore and aft limits 
of the lee side distribution and the keel side distribution are unchanged for 
these shifts (Fig. 4-2). 

The net pressure distribution for this load case is shown in Fig. 4-3 
for several stations along the length of the vehicle. This distribution consists 
of the variable external cavity collapse pressure and constant internal SRM 
ullage pressure on Branches 1 and 2 of the model. The loading on Branches 
3 and 4 consists of the variable internal and external cavity collapse pressure. 
Figures 4-3a, 4-3b and 4-3c show the net pressure distribution with the external 
cavity collapse pressure in the nominal position, shifted aft a distance equal 
to l/8 the SRM diameter and shifted aft a distance equal to l/4 the SRM diam- 
eter, respectively. All other values are constant for the three load distri- 
butions. The nozzle loads acting on the aft dome due to the hydrostatic 
pressure on the nozzle bell are small and are neglected in this final analysis. 

The expressions for determining the circumferential distribution of the 
internal and external pressures were coded for automatic generation of the 
pressures in Subroutine USRLD of STAGS. Listings of the coding for the 
nominal, l/8 shift and 1/4 shift distributions are given in Tables 4-1, 4-2 and 
4-3, respectively. 
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Fig.4-3a - Net Pressure Distribution 85/32/5, Nominal Position 
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Table. 4-1 

USRLiD LISTING FOR NOMINAL LOAD POSITION 
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310 t fA«=i.000Cb78olPL£C)»*Z.l + ^Cl'r63 I J»'P'L|.E" 
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Table 4-1 (Concluded) 


KC"(C0S(3.l<«1592*¥PHI I )**ETA 

PE-PKEEL— *Kt-«-<PtES-PKEEt7 

320 CONTINUE 

PIN- Pn»»S* tT»*C.0g tPH] 1T*TPTT"'PHTI — 

P-P l N- ( PE*FN ) 

GO”T(r70tr 

c ••• uKanCH h 

-H00 PL-RI.B 

Prt * I b « 5 

- p i n *apn +.5*u.“coytPHiiT*rPir-PH) 

PNEEL*33*6a6-,oHo*i«xL 

IRXL.6 E^35T3 rn"Gt>"'TtllTB 7 

PLtC»lb0i6l2-.H*»l»Xl, 

tO TO 920 

*H0 PLLC-201.6H-J .ea7*xu 

"' 1 *2o continue - ' 

IP (YOU .GT.90o) GO TO 130 

— P£« _ PKEEX — 

CO TO H90 

“ N 30 E T a « t 0 0 0 0 57 B * T p i. E ET« * 27 f+'t'OTS' JT •TCFE’ 
KC* (COS ( 3 t 1 ‘i 1 5926-PH I ) )*»ETA 

PE* PKEEL - + KX*TPL‘EE*PKtE:trJ 

<*40 P»PlN-(Pt«FN) 

900 CONTINUE — 

CALL fORCt(P,r»3iL t M> 

"99? CONTINUE 1 

RETURN 


SSf.SSf 
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Table 4-2 

USRLD LISTING FOR 1/8 LOAD SHIFT AFT 


sro?,’JrfS US9L0»USSLD 

zz -l a / 1 e./ 7 6 -o a F 58 .*o*-<-Q-rii- 

-2 

SJ390 JTINE -US7L0- 


c *** DATA F40H CM/ITT COLLAPSE V tf =8 5 , 3 / ¥H / 63 » -l 3/TH/1Q 
C***— LOAD -FOB— 85 /32^*! 

IFCt.ra.2J RETURN 
1EAL KC 

33 399 L = l,M'?oa 

33 999 4 = 1,>IC3L 

FH-1 . 1 7 S 

P^rrcrcm./lSO 

Xl~i.il J 

- - - IFC I3RSCH.E8.W GO- 1 

IF CI3RSCH.E3.21S0 TD 200 
... - IF { I39MCrl. E3.31GO -T0-30O- 

IF{I3RMCr).E8.4)S0 TO 400 

31 ANCH-1--*** — 

133 PIS=-t.S5 


. REV.D OAT* 


IFCXL.SE. 193. S> -60— TO— l 30- 
P<£El"23.+.3775*XL 
-• - IF (XL. SE. 56.51 GO— TO— 1 10“ 

PL EE =3. 3 

30 to 130 — 

113 IFCXL.GE.I35.5J SO TO 120 

- 9 LEE fi. 25 MXL- 56.51 

S3 TO 1 9 D 


120 a LlE ~103.-*-.964*CXL-l J6-*£l 

30 TO 190 

— 130 a XEEL = 4 2.9 25~.04O4-**rL 

IF CXL.3E.21l.75J 30 T3 140 

PL EE = 14 9.24 9*. 864*5 XL -193. 51- - 
33 ID 190 

143 PLEEfI 55. -.60 3* (XL -2 11. -7 54 

193 COStlSUE 

— IF(r(H1.5T.93 .3— GO — T-0 — 1 8 0 * - 

PETOKEELFFV 

33 10 173 .... 

150 r T AFj .3303578*1 PLEE1 **2.J*»Q1634*PLEE 

- XC -'CC3SC3. 1415926— P Hill* *ET* 

p E - 1 PX CEL «<CMPLEE-PXEEU >*FN 

— 173 -Pf PIH-OE 

so ro 900 

C ** 3RASCH 2 

230 PIMr-1.55 

• °L -4 5 . 5 

0‘C -1 5.5 
.._ — pf ; ( PL J 
P = 9IH -PE 
30 10 930 • 

C 9 RANCH 3 *** 

333 3 L F4 5 . 5 -■ -- 

?H :15.5 


F-XEELF34.12 4~.94-09^XL — — 

?LEEfj52«93-.s03*XL 
I f ( r ( >c 1 .ST. 90.1- SO TO 310- 


^IQUSIAT b 

°"«««2sr 
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Table 4-2 (Concluded) 


SO TO 320 

310 £TA=(. 3300576* tPLEE) *92. )*.01639*f»LEE' 

XC = tC0S(3.19irS»2fr-PHI7-»-**tT-*; 

• PE - PXEEL ♦XCMPLEE-PKEEU 

— S-JO-CONTISUE 

PI9= PH«-.5*tl.-C0$(PHt>)*{ PL -PH I ’ 

p=pIM-(PE*FV» 

30 TO 930 

C- 3 9 A NCH * ♦** 

935 3 L=91.5 

P H : 1 5 * 5 — — 

3 I9:P4 ♦.5*(1.-C0S(PHI11*<PL-PH» 

- PK£EL-33.bB6— s -3 904*X-L- 

IF(XL.GE.3b.3161 SO TO 910 

... - PLE£ = 155.3-.&33*XL- 

33 TO 923 

9-13 3 LEE-1 35^3-HH 

9 ? 3 CD'JTIHUE 

- - I F f r < *9 ) .3T. 90*1-60- TO-4-30 

3 £ - PKEEL 

• - • 50 TO 990 - -- 

933 LT A: f .3300S78MPLEEI **2. 1 ♦ .31634*PLEE 

9C -(COST-?-. 193-59 26 - PH I 1 > »*CT A 

PE= PKEEL ♦ XCMPLEE-PKEEL1 

993 3 1 V- < D E *FH l 

933 C0VTI9UE 

CALL. F09CE(P,9r3*t»-Hl- 

999 CtmiHUE 

2ETJ9S 1 
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Table 4-3 

USRLD LISTING FOR 1/4 LOAD SHIFT AFT 


k>f UK ,0'Ab USKLO.UShLD 

USA 1.3 »10/2J/> -tq, 1) 

-2 

o soL-hou r i wfc ubnuu 

c L>ATA FMJM CAVITY COLLAPSE V V« 0 S i 3 / VH /60 ■ - 1 0 / Th/ 10 HfV.O DATA 

i«*» 1.0 au H'« t»b/32/*b' CAVil Y-COLCAfSfc 'L-£rAfr“CTA«rTKA'*'' rREas»J M 

If IA.UQ. 2 J NLTUKN 

NC “ 

im vv? L*i»rmot. 

I Ilj V ? n*l»nc.uL * " 

F ^ *■ J V l 7 Cl 

PrlI = |T|Hl / Hi J . > * 3 • l'*! I 5V2« 

AL’ULl 

\Y i , oKHCH.uk. 1 iuU 10 100" ‘ — 

I r t i MS' NCH . ■- II > / ) u(J To 2 Utl 

it ( ii.AncM.bit . 1 1 oo Tu 3 Co ' — " 

I I 1 J :tkiiCn . Lu , •! I L ,0 To HOC 

L O K A t.Xn I »»« 

10 'J Kl».-I * 

l f t X l . 0 t ♦ 1 V .1 . b ) 0 0 To 11 0 — 

PnU!_l-23.t.o//0»XL 

II oO 10 1 2C ' ’ 

P L 1 1 -• J • O 

. 00 i o 1 vO — 

t J Ij Is l A L • o r. • J 3 6 . •} J 0 0 I i 1 1 3 tl 

PL - l. = 1 < 2b» t Xl -56 .’j J 

1 v 1*0 

I. *c* ll l.t.*’ I i'.l • .»oVsj^ l AL W 136 iS ) ' " 

GO fo 1 V 0 

"lit! Pil:.‘ l>t2.cii-j-.'.”iyR»At--- 

I’Ll t •- lJ 9.<,2 + . e , 7 b»lAU-lVJ,b) 

IF 1 r »n> .ol ,9S. 1 00 TO 180 

00 TO l-’O 

II. J C t • ^ 0 "J 05 /«» tPLut ) »» 2 . T+-S.J t 6 '^‘+»Pl.eE'-' 

i;o-=;(.uijfJ.l‘UbVi<>-t'iiill«»6TA 

PL-U’AutU " I PLt,t-PKKEU ) J »PN 

■| ? 3 P« Pl’l-Pt 

oU I -J VUO ' - 

L .» AiCs.iUit 2 *•« 

iutj Plrl = -1 .sS • 

Pi --IS . S 

PH 31 1 5 • 4 " 

Pu-U’u +<ij*iPrt-f't!*(i • “COS ( Ptl 1 ) ) i 
f‘t.Ptrt -Ft • - - 

0 Ll I IS V 0 3 

L ..p. L>>' A lit M •I-*..’ - * 

JoJ l*i.»Hb*b 

Pr:“ I j t j -- 

i , ..<-_l“J‘I. 12 S-, 0 H 0 H»AL 

PlI t ;| , J 1 '. 1 )*. 6 V 5 .*I. " 

1 r ( T ( it > • u T . Fu . ) 00 TO 31 Q 

l*L*t I.UtL ’• 

oo ro J2‘J 

»1U tT i M .OJUJt>; d* ( l‘LLU 1 **Z. ) ♦ * 0 liJ*t«PUeC’ 
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Table 4-3 (Concluded) 


K.C*(C0S(3»14l!i926«Pm J >»«ETa 

P£» HKEEU- *rc« TPLEt-PKEECr 

320 CONTINUE 

PIU" PH+0*U .-COi<PHtn«tPU-PH) 

P*P i N“ l Pt *K n ) 

fau TO <i 00'" 

C d si A N C rt 4 

luO Pl“H l .5 ■■ — - - 

f) H 4 I b » j 

PI '•■‘PH +»S*U .“(.OSlPHr) J*{PU-PHJ- 

P^Lfc.!.»33»4d6-.OHO‘<‘}lL 

■ I r < XL.QE . 35, 3 t'6 J ‘SO 'TTJ--420 

If I AL.C,d.3.dl6f :0 TO “10 
PL^L-162* J3+.6?b*Xl- 
faO [ O H VO 

4lo PLtE-> I o5*“ • f02« IXL-3.3 r* f 

uU TO 470 

420 i>£ELM3!b*-l.t337« ( XL-3ST3 f^T ‘ 

*t v u continue 

I F ( T ( ft j .GT. 7 3* ) ' 50 T0—4^a 

f u * i 

.,o ro 4 4'J -- 

43 J ETA-st.-jjJ0'j/j»{| , utt)-»*2,j + .0l63H*>Pl,EE 

* ICC * (COS ( J. I H l ->72 6 -PH I J ) «*ET.4' 

P»-’ PKUL «• K.C* tPLtL-PKEEU 

4 NO P- 1 ’ I >4- (f'C»Fi( I 

vu-j CO'i r i iiul 

CALI. F04CC(P,H.3»L,HI 

7 7 7 CJ’ifirfijL 

KcT'J iN “■ * * — — — — 

tNO 


0&f^4z b 
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4.2 SRB STAGS C LINEAR ANALYSIS 

Linear analyses were performed, for the 85/32/5 (VV/VH/9) case for 
both the original and rotated models. Each model was analyzed for all three 
pressure peak locations. No design factors were applied to the loads in the 
computer analyses. The. linear analyses showed good margins against the 
material yield allowables in all structural areas of the SRB. Figures 
4-4 and 4-5 show the maximum stress values determined in the linear 
analysis for the three pressure locations. The first stress value is for the 
max pressure peak in the nominal position, the values denoted (1/8) and (1/4) 
are for the pressure peak shifted aft 1/8 D and 1/4 D, respectively. 

From Figs. 4-4a and 4-5a the maximum stress values in the SRM case 
and dome structure are seen to be well below the material allowables for the 
DbAC steel. If a design factor (DF) of 1.25 is applied to the maximum stress 
value given for the dome skin, f = 1.25 (14.6, 400) = 183,000 psi, the resultant 
stress value still falls below the ultimate material allowable (Table 3-5) giving 
a positive margin of safety. 

The maximum stress in the steel reinforcing tee -ring at the cylinder - 
conical skirt interface is seen to be well below the maximum allowable stress 
value for both model orientations, Figs.4-4b and 4-5b. 


The computed stress values in the skirt stucture for the linear analyses 

% 

are also shown in Figs.4-4b and 4-5b. These maximum values are seen to 
be well below the material allowables for the 2219 aluminum. The largest 
stress value occurs in the skin just forward of the upper ring at (f> = 180 deg, 
lee side, for original model configuration. Applying a design factor of 1.25 
produces a stress value still below the ultimate material allowable. Hence, 
the margin of safety is 


MS = 


tu 


(DF)(f) 


- 1 


63,000 

(1.25)(46,200) 


1 = + .09 
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Stress in psi 
r = ring, s = skin 
c = compression 
f * tension 



Fig.4-4a - Maximum Stresses, Original Model, 85/32/5 Linear Analysis 
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Note: Heavy lines denote stiffeners. 


Fig. 4-4b - Maximum Str 
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Stress in psi 


NOTE: Heavy lines denote 
stiffener &. 


t = ring, s = skin 
c = compression 
t = tension 



Eig.4-5a - Maximum Stresses, Rotated Model (+_Zg) 85/32/5 
Linear Analysis 
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Sta. 1828.850 

Branch 3 
Sta. 1839.684 


Branch 4 


I 

J— 4 
vO 



f = 51,310 

50,890 (1/8) 
50,930 (l/4>: 


\ Upper Ring 
;■ f = 20,160 I 

tr 20,180 (1/8)! 
20,290 (1/4) 
f = 35,760 


35,300 (1/8)' 
35,070 (1/4): 


. f tr = 17,260 

17,410 (1/8); 


17,640 (1/4)1 


Fig.4-5b - Maximum Stresses, Rotated Model (+Zg) 85/32/5 Linear Analysis 
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The maximum stresses in the rings are also given in the figures. These 
values should be checked against the flange crippling stresses in a detail 
stress analysis to ensure structural integrity. From the size and cross- 
sectional dimensions of the rings, it would appear by inspection that local 
crippling would not be critical but this should be determined by calculation. 
The displacements normal to the surface for the $= 0, 90 and 180 deg meri- 
dians are shown in Figs. 4-6 and 4-8. Figure 4-7 shows the deformed shape 
at several cross-sectional locations for the original model, the rotated model 
has a similar general deformation pattern. 

4.3 BIFURCATION BUCKLING ANALYSIS 

Bifurcation buckling analyses were performed for the 85/32/5 load case 
for both the original and rotated models. The results of these studies are 
given in Table 4-4. The eigenvalues are seen to be at least 2.5 times the 
applied load with the predicted buckle occurring in the cylinder at station 
1797.8 between the aft tee ring and clevis joint for both models. The buckling 
mode shape is shown in Fig. 4-9. Even with an appropriate knockdown factor 
■ (KDF) of 0.63 applied (Ref. 3), the resultant predicted buckling load factor (LF) 
is well above the design load. 

LF - KDF x Eigenvalue 
LF = (0.63) (2.526) = 1.592 

This factor results in a. positive margin of safety even with the 1.25 design 
factor applied* to the loading. 


MS = 



1.592 

1.25 


1 = + 0.274 


4.4 NONLINEAR ANALYSIS 

A nonlinear analysis was performed to more accurately determine the 
failure load. The results of this analysis at Station 1797.8, 4 > ~ 180 deg is 
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Fig. 4-6 - 85/32/5 Original Model (+ Yg) Deformed Shape, Nominal 
Load Position 


4-21 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


Normal Displacement {In.) 


LMSC-HREC TR D497054 



Fig. 4-7 - Normal Displacements, Original Model, 
85/32/5 Nominal 'Position 
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Fig. 4-8 - 85/32/5 Rotated Model (+Z R ) Deformed Shape, l/8 Diameter 
Aft Load Shift 
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Table 4-4 

BIFURCATION BUCKLING, 85/32/5 LOAD CASE 


SRB Model/Load 
Po sition 

Eigenvalue 

Location 

Original Model (+Yg) 

' 


Nominal 

2.569 

Sta. 1797.8 @ 180 deg 

1/8 Aft Shift 

2.785 

Sta. 1797.8 @ 180 deg 

1/4 Aft Shift 

2.993 

Sta. 1797.8 @ 180 deg 

Rotated Model (+Z_) 
— ’ ±3 



Nominal 

2.526 

Sta. 1797.8 @ 180 deg 

1/8 Aft Shift 

2.750 

Sta. 1797.8 @ 180 deg 

1/4 Aft Shirt 

2.968 

Sta. 1797.8 @ 180 deg 
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shown in Fig. 4- 10 for the rotated model with the 85/32/5 load case in the 
nominal position. The linear, nonlinear and bifurcation buckling analyses 
are compared in this figure. The nonlinear collapse is seen to occur at a 
load factor of approximately 1.8. This factor is seen to be well above 
the applied load, LF = 1.0, indicating structural integrity will be maintained. 


Stress values in the skirt structure for the applied load (load factor of 
1.0) are shown in Fig. 4-11 for the nonlinear analysis of the rotated model. 
The major skin stress value near the Upper Ring at $ = 180 deg is seen to be 
25% greater than the corresponding linear analysis value given. in Fig. 4-5b 

{ = 1.25). The stress values in Fig. 4-11 for the rotated model skirt 

are seen to be below the material yield allowable. With the design factor of 
1.25 applied, the maximum stress value is still below .the material ultimate, 
f = 44,800 (1.25.) = 56,000 < F t . However/ if this ratio of nonlinear/linear 
value of 1.25 is applied to the maximum stress value (f = 46,200 psi) for the 
original model linear analysis shown in Fig.4-4b, the maximum skin stress 
value will be f = 57,750 psi. Using the design factor of 1.25 the resultant 
stress is greater than the material ultimate giving a negative margin of 
safety in the skin at this area of the skirt. 


MS 


tu 


(DF) f 


1 


63,000 

" (1.25)(57,750) 


MS = -.127 


A more detailed analysis should be performed for this area. Other areas of 
the structure have positive margins of safety by inspection. 
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Fig. 4-10 - SRB Rotated Model, Sta. 1797.8, Lee Side, 85/32/5 
Nominal Position 
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Note- Heavy lines denote stiffener b. 


Fig. 4-11 - Maximum Stresses Rotated Model, 85/32/5 Nonlinear Analysis 
Nominal Load Position 
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4.5 THERMAL, EFFECTS 

To determine the sensitivity of the aft SRB structure to thermal loads, 
the rotated model described in Section 3 was used. However, the discrete 
stiffener used in previous models could not be used. To account for thermal 
gradients each stiffener cross section must be divided into sub-elements so 
the thermal expansion may be considered constant within each sub-element. 

Stiffeners located on the exterior of the motor case were divided into 
4 or 5 sub-elements so that temperatures could be allowed to vary linearly 
from the outer portion of the stiffener to the skin temperature at the sub- 
element nearest the skin. For internal stiffeners (those on the skirt or at 
the aft end of the dome) only one subelement was used since the stiffener and 
adjacent skin were of the same temperature. No variations in temperature 
through the skin thickness were allowed. 

The skin and stiffener temperatures given in Table 4-5 were used in 
the analysis. A stress-free temperature of 70 F was used. 


Table 4-5 

TEMPERATURE DISTRIBUTIONS FOR THERMAL 
EFFECTS DETERMINATION 



Motor Case 

Dome 

Skirt 


Temperature (F) 

Skin 

250 

250 

' 160 

Stiffener 

400 

250 

160 


Maximum stresses resulting from the 1/4 D shift pressure distribu- 
tion described previously combined with the temperatures above are shown 
in Fig. 4-12. The skin stresses in the skirt at $ = 180 deg, adjacent 
to the upper ring is the only area within the skirt where problems may arise. 
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i 

w 

o 


Sta. 1828.850 
Branch 3 
Sta. 1839*684 


Branch 4 


Cone-CYL Intersection, 
Reinforcing Ring 


Upper Ring 

f = 20,070 
cr 

f ts = 39,810 



Intermediate 

Ring 

£, . = 19,500 


Aft Ring 


Note: Heavy lines denote stiffeners 


Fig. 4- 12 - Maximum Stress Rotated Model, 85/32/5 Plus Thermal Distribution 
1/4 Aft Shift 
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The stress in the skin increased by approximately 4740 psi due to the inclusion 
of thermal’ effects. For the rotated model this increase does not indicate diffi- 
culties. However, the original model could be expected to increase by the same 
amount. In Fig.4-4b it was shown that the maximum skin stress was 46,200 
psi for the original model. An increase of 4,740 psi increases the maximum 
stress to 50,940 psi. Applying a design factor of 1.25 to this value the stress 
becomes 63,675 psi, approximately the ultimate stress for the 2219 aluminum 
used in the skirt. 
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5. REVIEW OF BROWN ENGINEERING SRB STAGS MODELS 

Brown Engineering personnel, under support contract to NASA-MSFC, 
experienced several difficulties in obtaining accurate results using STAGS. 
These problem, areas were investigated and causes were determined. Most 
problems were associated with the grid network chosen and the ordering of 
the branches of the model. The particular modeling problems were covered 
by the parameter studies and other grid studies in Section 2. 

In the Brown SRB STAGS math models that were reviewed, the cone 
branch is designated as Branch 1 and the cylinder as Branch 2. The cone 
longitudinal grid originates at the intersection and terminates at the aft end 
of the skirt. The cylinder grid originates at the forward end of the cylinder 
segment (support plane) and terminates at the cone-cylinder intersection. 
This is an inefficient method for the ordering of the terms in the handed 
matrix solution routine used in STAGS. The most efficient method is for 
the cylinder to be designated Branch 1, originating at the forward end of the 
cylinder segment and terminating at the cone -cylinder intersection. The 
conical skirt would be Branch 2 with the origin at the cone cylinder inter- 
section and terminating at the aft end of the skirt. This provides a con- 
tinuous flow of node connections throughout the model, thus reducing the 
number of interconnected matrix terms. (This is the ordering sequence 
used in the model constructed in Section 3.) 

A large savings in computer time is realized using this orientation. 
Table 5-1 shows a comparison of Univac 1108 computer run times for com- 
parable Brown and Lockheed STAGS B SRB models with a cavity collapse 
load case. The ordering sequence shown in Section 3 was used for the four- 
branch Lockheed model. A saving of 2 hours and 6 minutes of charged time 
is realized by the proper ordering sequence. 
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Table 5-1 

SRB STAGS B MODEL COMPARISON ' 


Item 

Lockheed 

Model 

Brown 

Model 

Number of Branches 

4 

2 

Number of Nodes 

3183 

2394 

Number of Equations 

5187 

5786 

CPU Time* 

0/19/21.689 

0/36/46.987 

I/O Time 1 

0/6/38.448 

1/53/26.078 

Total Time V 

0/27/2.868 

2/32/52.909 


^Time is hours/minutes/seconds. 


The Brown models have 21 finely spaced grid divisions in the longitudinal 
direction on the skirt but only eight sparsely spaced divisions on the cylinder. 
From Fig. 2-6 of Section 2 this is shown to inhibit the true deformations at the 
cone-cylinder intersection and results in higher predicted stress- values in the 
reinforcing ring. This grid spacing would also result in higher computed values 
for the buckling load. The accuracy of these models could be improved by add- 
ing several closely spaced grid lines adjacent to the intersection and reinforcing 
ring on the cylinder branch. 

The SRM aft dome cap has not been included in the Brown models. The 
addition of this structure would give radial stiffness to the cylinder wall in the 
area of the highest pressure in the cavity collapse load. Also the constant 
internal pressure load on the skirt, dome and cylinder was not included in 
the load data. These corrections should be made to obtain accurate results. 
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6. CONCLUSIONS AND RECOMMENDATIONS 
6.1 CONCLUSIONS 

The STAGS C finite difference computer program is an efficient and 
accurate tool for the structural analysis of general shell structures. Param- 
eter studies should be performed and several grid patterns investigated using 
a simplified model of the structure and loading before the large detailed model 
is constructed. This method will produce a mathematical model that has the 
most efficient computer run time and one that will give accurate results for 
the linear, nonlinear and bifurcation buckling analyses. 

The 86/32/5 (w/VH/ 0) cavity collapse load case was determined to 
be the maximum loading for the SRB skirt structure and the aft portion of the 
motor case. The linear analyses showed that structural integrity will be main- 
tained for the applied load. With a design factor of 1.25 the linear analyses 
stresses will still give positive margins of safety using the material ultimate 
allowable. The bifurcation buckling analyses predicted buckling in the cylinder 
motor case at Station 1797.8 on the lee side of the structure. This location 
was typical for both the rotated (+Zg) and original (+ Y^) models for all three 
load positions. The minimum eigenvalue calculated was 2.526. Even with an 
appropriate knockdown factor of 0.63 applied, the resultant eigenvalue of 1.592 
gives a positive margin of safety for the structure. 

The nonlinear analysis at the applied load (LF = 1.0) show stress values 
in the inelastic range in the skirt for the original model (+Yg). This occurs 
in the skin just forward of the upper ring at the lee side. Applying the design 
factor of 1.25 results in a stress value above the material ultimate giving a 
negative margin of safety {- .127) at this location. Detailed analyses should 
be made of this area to more accurately predict the true inelastic stress 
distribution and the corresponding margin of safety. 
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6.2 AREAS FOR ADDITIONAL STUDY 

Some areas for future study that need to be investigated and incorpo- 
rated into the STAGS SRB mathematical model are the addition of the pressure 
normal to the web of the skirt rings. This will have negligible effect on the 
overall skirt structure but may significantly affect the stress levels in the 
ring elements. For the thermal analysis of the skirt, a constant temperature 
was used. A more accurate temperature distribution over the cylinder, dome 
and skirt as well as thermal gradients through the skin of the structure should 
be obtained and a combined pressure/thermal analysis made. Structural anal- 
yses should be performed for the cylinder for the pressure peak shifted forward 
1/8 and 1/4 SRM diameters. Linear, bifurcation and nonlinear analyses with 
thermal effects, included should be performed. A simplified 360 deg model 
should be coded and analyzed with the lee and keel meridians off the planes 
of symmetry. These results could be compared to the analysis of the same 
model with the loading on the symmetry planes. A comparison of these 
models would give a factor to be applied to the 180 deg detail model stresses 
to obtain relative stress magnitudes. The addition of a plot package to check 
data input and structure geometry would be a very beneficial addition to the 
STAGS C program. The deformed structure plots would add greatly to the 
efficient utilization of this program during structural analysis tasks. 
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